Journal of Theoretical Biology

Reputation plays a major role in supporting cooperation among unrelated individuals through indirect reciprocity. By helping others, individuals build a good personal reputation and receive greater benefits from future partners. Most models of indirect reciprocity assume that a person's reputation reflects only their own behaviour. Yet in many societies, people are also judged by their family's reputation. How family reputation affects the evolution of cooperation, and whether reliance on it can itself evolve, remain unclear. Here we show that reputation inheritance expands the conditions under which indirect reciprocity favours cooperation, increasing helping and favouring greater reciprocity. Greater reciprocity in turn favours stronger reliance on inherited reputation, creating a positive feedback that stabilises cooperation, especially when interactions are infrequent or personal behaviour is difficult to observe. This feedback arises because cooperation generates future benefits both for the individual, through their personal reputation, and for their descendants, through inherited reputation. Reputation inheritance thereby provides a route via which kin selection and reciprocity, often treated as alternative explanations for cooperation, can reinforce one another. Our model helps explain why family-based reputation occurs across diverse human societies and provides an evolutionary framework for studying phenomena organised around family standing, including kin-based institutions, feuds between families and honour-based violence within them.

Inappropriate antibiotic use presents a major global health challenge, particularly in low-resource settings where access to quality care is limited but antibiotics remain relatively unrestricted. This study estimates the causal effect of frontline primary care quality on inappropriate community antibiotic use, combining detailed community-based data from approximately 100 rural villages in rural China with an instrumental variable (IV) approach embedded within a double/debiased machine learning (DML) framework. We linked objective measures of village doctor clinical practice quality, measured through unannounced standardized patient visits, to household-level antibiotic use data collected from the same villages. To identify the causal effect, we constructed multiple candidate instruments from extensive provider characteristics and used an ensemble of machine learning algorithms within a flexible DML-IV framework to approximate an optimal instrument, addressing a many-weak-instruments problem. We found that improving village provider clinical practice quality reduced both antibiotic receipt during healthcare encounters for common diseases and household antibiotic storage for future self-medication. Our findings suggest that strengthening frontline primary care quality can meaningfully reduce inappropriate community antibiotic use without restricting access to essential treatment. More broadly, this study illustrates how causal machine learning can strengthen conventional causal estimation in complex observational settings in global health economics research.

Introduction Hearing loss (HL) affects over 1{middle dot}5 billion people globally and India shares a disproportionately high burden including Disabling Hearing Loss (DHL). HL affects an Individual socio-economically, but there are limited studies on the broader societal economic consequences of HL in India.Methods Using Cost-of-Illness (COI) approach, we studied the societal economic burden of HL in India. This study uses epidemiological and macroeconomic data and modelling to estimate the loss of Gross National Income (GNI) due to HL and DHL across three economic pathways. Uncertainty is evaluated using deterministic and Probabilistic Sensitivity Analyses (PSA).Results The model estimates that there are in India, 289 million and 85{middle dot}9 million people with HL and DHL respectively. Direct Loss of GNI and Indirect Loss of GNI (Caregiver burden) are estimated as INR 4,648{middle dot}4 billion (USD 55{middle dot}6 billion) and INR 3,268 billion (USD 39 billion) respectively. The Loss of GNI due to Low Education amongst those with HL is estimated as INR 1,041{middle dot}9 billion (USD 12{middle dot}45 billion).Discussion Economic burden of HL is presented across three pathways with Direct Loss of GNI due to DHL being the greatest. It also presents age stratified caregiver economic burden. The findings of the study help in estimating similar cost pathways, advocacy, and policy decisions towards reducing HL prevalence in India and LMICs. This study also highlights the need for India specific estimations related to the HL attributable low education, state-wise disaggregates, and prevalence studies. Funding This study has not received any funding.

Background Understanding the mechanisms underlying brain aging and age-related pathological changes is essential for advancing brain health research. Our group previously developed a mechanistic mathematical model of healthy brain, Chamberland et al. (2024) that integrates key biological processes involved in normal aging, from which Alzheimer's disease (AD) related changes may emerge naturally. Objectives To characterize and validate this brain model by evaluating its sensitivity, calibrating its parameters, and assessing generalizability in independent populations. Methods The model represents the evolution of key biological processes associated with brain aging, including amyloid beta (A{beta}), tau pathologies, neuroinflammation, and neuronal death. After identifying the 30 most influential parameters, we calibrated the model using cognitively normal (CN) participants from the AD Neuroimaging Initiative (ADNI) database (n = 211) by minimizing a loss function composed of three outcomes (AB) plaques, tau tangles, and neuronal density). The calibrated model was then applied to the UK Biobank cohort (n = 35,899) of normal controls (aged 44-82 years). The effects of sex and APOE were evaluated using stratified simulations. Results Parameter calibration significantly reduced the prediction errors for A{beta} and tau. Neuronal density predictions showed strong agreement in the UK Biobank cohort. The variance decomposition identified APOE status as a major contributor to variability in A{beta}. Conclusion Our validated brain health model links mechanistic pathways with population data and reproduces neuronal density patterns in an independent cohort. These findings support its use as a framework for studying brain aging and investigating how Alzheimer's disease related pathological changes may emerge with aging.

Alzheimer's disease is characterized by the accumulation and aggregation of amyloid-{beta}(A{beta}), but the molecular mechanisms linking environmental and infectious factors to A$\beta$ conformational changes remain incompletely understood. Herpes simplex virus type 1 (HSV-1) has been proposed as a potential contributor to AD pathology, and interactions between the viral glycoprotein B (gB) and A$\beta$ may influence the conformational behaviour of the peptide. Molecular dynamics (MD) simulations provide atomic-scale information on such interactions, but conventional structural descriptors may not fully capture changes in the organization of residue interaction networks. Here, we introduce a graph-geometric framework based on Forman-Ricci curvature to characterize the evolution of residue interaction networks during MD simulations. Each simulation frame is represented as a residue interaction graph based on C--C contacts, and residue-wise curvature profiles are analysed across time. We apply the framework to A{beta}1-42 in isolation and in complex with HSV-1 gB. Conventional MD analyses indicate stable association of the simulated complex, favourable interaction energetics, and conformational changes in A{beta}, including a transition from -helical structure toward {beta}-turn-rich conformations over the simulated timescale. Forman-Ricci curvature reveals pronounced and spatially localized remodelling of the A{beta} residue interaction network in the complex, with the strongest changes concentrated in the C-terminal region. These regions also exhibit reduced temporal curvature fluctuations and progressively distinct geometric behaviour throughout the simulation. Hierarchical clustering further identifies cooperative groups of residues with coordinated curvature dynamics, including a prominent C-terminal domain. Together, these results demonstrate that Forman-Ricci curvature provides a complementary description of biomolecular dynamics by capturing changes in the geometric organization of residue interaction networks that are not directly represented by conventional structural descriptors. The framework provides a general computational approach for studying network-level structural remodelling in protein molecular dynamics and offers a quantitative perspective on the conformational consequences of HSV-1 gB--A{beta} association.

Vaccine effectiveness (VE) estimates can vary widely between years and populations, even for the same vaccine. Estimated VE is known to be sensitive to susceptible depletion and differences in pre-vaccination infection risk between vaccinated and unvaccinated populations. However, how variation in pre-vaccination risk within and between the two groups affects VE estimates over time remains unclear. This uncertainty is especially important given negative VE estimates. We investigated the difference between estimated VE and true vaccine protection considering continuous distributions of pre-vaccination infection risk under three scenarios. When the vaccinated and unvaccinated populations differ in their mean risk, estimated VE can be higher or lower than true vaccine protection. Similar patterns arise when both populations share identical means but different risk distributions. Finally, if infection-derived immunity lasts longer than vaccine protection, annual VE estimates can vary by tens of percentage points between years despite constant true vaccine protection. These theoretical results underscore that VE studies estimate contrasting risk between vaccinated and unvaccinated individuals in a particular time and place, and VE estimates can vary counterintuitively between years and populations even with constant vaccine-induced protection. Explaining variability in estimated VE thus requires a more complete understanding of populations' distributions of infection risk.

A dental AI system and a dentist reading the same radiographs form a paired comparison. Published comparative studies often report the two arms separately against a reference standard, leaving the joint pattern of correctness between them unavailable for secondary paired inference. We show what that omission costs. The accuracy difference remains exactly identified; its sampling variance does not, so the report contains the estimate and not its uncertainty. On a study of 282 units, two published accuracies are consistent with 38 distinct joint tables whose confidence intervals differ in width by a factor of 2.5. The consequence is a three-zone decision map rather than a single threshold: differences at or below 1.06 points are non-significant under every compatible table, differences at or above 6.03 points are significant under every compatible table, and in between the published numbers cannot decide. We then show the omission is repairable at negligible cost. One additional integer, the number of units both arms classify correctly, identifies the joint table exactly and restores standard paired inference. For a panel of readers the pairwise dependences must arise from one joint distribution, a constraint that binds once three readers are present; publishing each reader's joint-correct count against a single reference reader cannot widen and may tighten every pairwise bound, and in a 7-arm experiment reduced them by a median of 37% even for pairs excluding that reference. Where the integer was never published we give DentalPair-Cert, an interval with finite-sample coverage uniformly over every admissible within-unit AI-dentist dependence under the independent-sampling-unit model, certified in both the nuisance maximization and the inversion. Across 4,200,000 simulated comparisons an independence analysis falls to 74.5% coverage with 12.2% type-I error; in a purposive sample of 9 recent comparative studies, 1 reported a paired test on discordant units.

Contact inhibition of proliferation (CIP) enables epithelial tissues to self-regulate growth and maintain tissue homeostasis. However, how cell-level mechanical contact, tissue-scale structural order, and proliferation kinetics interplay remains a fundamental open question in living matter physics. Here, we present a particle-based model of a confluent epithelial monolayer governed by overdamped dynamics, where individual cells interact via a two-dimensional hard core- soft shoulder potential. By comparing structural evolution during quasistatic densification with previously reported experimental division kinetics, we find that the dynamics of proliferation arrest mimics the onset of direct steric contacts between the hard cores of the shell. Identifying hard core contacts as the physical driver of CIP, we couple our mechanical model with a stochastic Monte Carlo division scheme in which the instantaneous division rate decreases to zero from an intrinsic value as the number of hard core contact increases to six from zero. We demonstrate that for high intrinsic division rates, the cellular densification outpaces mechanical relaxation. This kinetic mismatch drives premature hard-core contact formation, shifts the onset of jamming and contact inhibition to lower packing fractions, and induces increasingly disordered transient configurations before the tissue universally converges to a hexagonal close-packed limit. Our model's predicted division kinetics and structural order evolution are consistent with epithelial monolayer experiments, both reported and our own. This minimal physical framework links single-cell steric contact mechanics directly to tissue-scale growth regulation and structural evolution.

Background: Lumbar disc herniation, with or without concomitant disc degeneration, is a major cause of lumbar radiculopathy and low back pain, which also a key public musculoskeletal disorder without an exact pathophysiology. Studies have suggested that inflammatory cells and biochemical markers of inflammation also play an important role in lumbar radiculopathy in addition to nerve compression. The aim of the present study was to assess the association of selected circulatory inflammatory markers (CRP, hs-CRP and E-selectin) in patients with lumbar disc herniation without radiological degeneration (LDH) and lumbar disc herniation with radiological degeneration (LDHD). Materials & methods: This case-control study included 208 participants, comprising 104 patients with lumbar disc pathology and 104 controls. Patients were further stratified into LDH (n=67) and LDHD (n=37). Serum CRP, hs-CRP and E-selectin concentrations were measured. Results: Among the patients, 35.6 % presented with LDHD while 64.4 % had only LDH. Significantly increased median hs-CRP (p<0.001) and CRP (p<0.001) were observed in patients groups compared to controls, while CRP showing a consistent independent association across the combined disease (OR=1.68, 95% CI=1.33-2.14, p<0.001), LDHD (OR=1.62, 95% CI=1.16-2.20, p=0.005) and LDH (OR=1.69, 95% CI=1.30-2.20, p<0.001) multivariable models. No significant difference was observed in serum E-selectin between the study groups. Multivariable models incorporating inflammatory and clinical variables demonstrated substantially greater discriminatory performance than individual biomarkers alone. Conclusion: Elevated circulating CRP and hs-CRP concentrations were associated with lumbar disc pathology, with CRP showing a consistent independent association across the combined disease, LDH and LDHD multivariable models, whereas E-selectin showed no significant association. Multivariable models incorporating inflammatory and clinical variables demonstrated greater discriminatory performance than individual biomarkers. These findings support a potential systemic inflammatory component in lumbar disc pathology, although the cross-sectional nature of the measurements does not establish causality or a local inflammatory response within the disc.

Aptamers are structured nucleic acid ligands capable of high affinity, high specificity molecular recognition generated using variations of the SELEX (Systematic Evolution of Ligands by Exponential Enrichment) process. However, SELEX often produces sequences that enrich yet may lack binding efficacy. We propose a measure called the Ruggedness Composite Index (RCI) along with a method for computing it, that can be used to distinguish binding-competent ('active') aptamers from weak or non-binding ('inactive') aptamers. Given a set of aptamers, RCI incorporates information on their fragmentation (landscape partitioning), basin entropy (metastable state distribution), cumulative density irregularity (non-uniform occupancy), and structural energy correlation length (structure-energy coupling scale). We test whether secondary-structure folding energy landscape topology distinguishes active from inactive aptamers using a multiscale level set framework across six datasets. Active aptamers show lower RCI values and occupy smoother, funnel-like conformational spaces, while inactive aptamers show higher RCI values, reflecting fragmented, high-entropy landscapes. By contrast, classical thermodynamic features, such as minimum free energy, show limited discrimination between active and inactive aptamers. In all datasets, sequences that exhibit enrichment which is not monotonic but lack specificity exhibit elevated ruggedness, indicating landscape topology can predict non-specific enrichment. These results indicate that folding landscape organization can be used as a predictor of aptamer activity and establish RCI as a simple, mechanistically interpretable measure for improving candidate prioritization, especially in therapeutic aptamer discovery.

Standard polygenic risk scores (PRSs) are constructed based on additive genome-wide association study (GWAS) summary statistics. Nonlinear machine learning methods have been increasingly applied to construct PRSs directly from individual-level data, with the aim of improving predictive performance over standard PRSs through their ability to model non-additive genetic effects. However, their superiority across studies has been inconsistent, and the conditions under which they provide meaningful improvements remain unclear. We combined theoretical analysis, simulations and a real-world application to investigate when two widely used nonlinear machine learning methods, random forest and XGBoost, outperform standard PRSs. Theoretical analysis showed that standard PRSs can implicitly capture part of the genetic variance attributable to nonadditive genetic effects through their contributions to marginal SNP effects, thereby losing less information than commonly assumed. Although nonlinear models have a higher theoretical potential, their greater flexibility incurs a bias-variance trade-off that can limit predictive gains at finite sample sizes. Simulations showed that XGBoost outperformed the standard PRS only when the genetic architecture involves a sufficiently large proportion of interaction genetic variance concentrated across relatively few interaction effects and large training samples were available. Random forest consistently underperformed the standard PRS. In an application to ischemic heart disease prediction using UK Biobank data, XGBoost showed no meaningful improvement in predictive performance over the standard PRS, whereas random forest again performed worse. Together, these findings suggest that nonlinear machine learning do not uniformly outperform standard PRSs; rather, their relative performance depends jointly on genetic architecture and training sample size. Our study helps to reconcile the inconsistent results reported across previous studies and provides a framework for identifying settings in which more complex PRS models are likely to be beneficial.

The erythrocyte sedimentation rate (ESR) is one of the most common and widely used laboratory diagnostic parameters in connection with inflammatory reactions and it is probable that every reader has already experienced a determination of their ESR. A rapid ESR is a non-specific parameter that provides information about the inflammatory process. Although the origins of this methodology date back to antiquity, the description of the process as the collapse of a percolating gel formed from erythrocytes has only recently been achieved. It was not yet known whether slow ESR has any medically relevant significance. Here we show a variety of clinical pictures that exhibit a systematically slow ESR (e.g., sickle cell disease, neuroacanthocytosis syndromes, chronic mountain sickness). Using a combination of measured data and physical modelling, we show how the accuracy and significance of ESR data can be increased. With this improved ESR (supraESR), we introduce a completely new, cost-effective diagnostic parameter, based on an established and easily automated measurement method, that enables low-cost screening for neuroacanthocytosis syndrome, a group of rare neurodegenerative diseases previously detectable only through complex diagnostic tests.

Background: Pulmonary arterial hypertension (PAH) is a serious disease with poor prognosis, especially in infants or preterm babies and there is still no optimal treatment for this disease. Noradrenalin (NE) is a vasoactive mediator which is released by sympathetic ganglion. According to previous studies, NE/1-adrenoreceptors is not only in regulating normal physiologic responses, but also in the pathogenesis of PAH. However, the mechanisms of NE in PAH are not fully understood. Methods: Human PASMC (PASMC) was used in this study. Cell viability assay and Wound healing assay were used to evaluate the proliferation and migration of PASMC. Immunoprecipitation and western blots analysis were used to investigate the mechanisms which involved in NE-induced PASMC proliferation. Results: We investigated that NE could induce human PASMC proliferation and migration. Furthermore, we first find that endothelin 1 (ET-1) signaling pathway plays an important role in NE-induced PASMC proliferation. ET1 is a critical molecular which is known for regulating cell growth and migration. We investigated that NE could increase NE-1 secretion, further enhancing ET-1 bind to its receptors. For further clarifying the downstream signals in NE/ET-1 induced PASMC proliferation, we detected the phosphorylation and expression levels of ERK and JNK. Conclusions: By combining the results from ours and previous studies, we believed that JNK/c-jun pathway may play an important role in NE-induced PASMC proliferation. Key Words: Noradrenaline; Pulmonary Arterial Hypertension; Pulmonary Artery Smooth Muscle Cells; Endothelin-1; JNK/c-Jun Signaling.

Bayesian borrowing of external information can improve trial efficiency, particularly in pediatric and rare disease settings where patient populations are limited, but may introduce bias and inflate the Type~I error rate when the trial differs from external studies. Recent U.S. Food and Drug Administration (FDA) draft Bayesian guidance emphasizes careful evaluation of external information, prior specification, and assessment of operating characteristics. This paper compares three meta-analytic-predictive (MAP)-based methods for Bayesian borrowing: the MAP prior, robust MAP (RMAP) prior, and self-adapting mixture (SAM) prior. An adaptive platform trial design in Rett syndrome is used as a case study. Simulation studies evaluate frequentist operating characteristics under varying prior--data conflict, between-study heterogeneity, treatment effects, and clinically significant differences (CSDs) for the SAM prior. The MAP prior achieved the greatest efficiency when external and current data were compatible but exhibited the largest bias under substantial prior--data conflict. The RMAP priors improved robustness through fixed robust-component weights, whereas the SAM prior adaptively adjusted borrowing and was less sensitive to prior--data conflict while retaining efficiency gains when the data were compatible. Although the CSD influenced the degree of adaptive borrowing, as reflected by effective sample size, it had only a modest impact on frequentist operating characteristics. Sensitivity analyses using a skeptical robust component yielded similar qualitative conclusions, while accentuating the differences between the MAP and RMAP priors. These findings provide guidance for evaluating and selecting MAP-based borrowing strategies before trial implementation, particularly in rare disease settings, consistent with current FDA recommendations.

The study of evolution between species (phylogenetics) and the study of evolution within a species (population genetics) are highly related, as the same biological mechanisms are fundamental to both fields. Although both have been studied for a long time, their joint study in a unified setting has been prevented by the different time scales they consider and the different data types they employ. A similar discrepancy holds for their whole-genome specializations, comparative genomics and pangenomics. Two active areas in these fields are genome rearrangement studies and graphical pangenomics, respectively. Since the emergence of graphical pangenomics, these have existed as separate fields, despite observations that central data structures representing genomic variants in both fields are highly similar. While there exists a wealth of theoretical results for various rearrangement models in comparative genomics, the application to pangenomic data is hampered by the limitations of rearrangement problem formulations. On the practical side, pangenomes typically contain too many individual genomes for classical problems, such as the often NP-hard parsimony problems, to be solved, or for all-vs-all comparisons using rearrangement distances to be performed. On the theoretical side, some assumptions in the formulation of rearrangement problems, such as the assumption of an underlying tree, are inadequate for many pangenomes. In this work, we propose the Complete Ancestral Reconstruction for Pangenomes (CARP) problem, which overcomes these limitations while retaining intuitive relationships to both classical rearrangement problems and pangenome graphs.

This study aimed to quantify the longitudinal tracking and cross-sectional construct validity of a single-item questionnaire measuring recreational physical activity frequency (RPAF) in the Womens Healthy Ageing Project. At baseline, 474 participants aged 45-55 reported RPAF from 1993 to 2014. Longitudinal tracking of the RPAF item was assessed as a consecutive-wave and baseline-referenced measure using linear weighted kappa (LWK), Spearman correlations, exact agreement and within-one-category agreement. Construct validity in the form of convergent and known-group validity was assessed using the International Physical Activity Questionnaire (IPAQ) leisure activity domains, Short Form 36 physical function (SF-36-PF) subscale, Timed Up and Go (TUG), hand grip strength (HGS) and waist-to-height ratio (WHtR). 474 participants provided baseline RPAF data. Pairwise longitudinal samples ranged from 176 to 459 across the study. Consecutive-wave LWK ranged from 0.38 to 0.49, and Spearman correlations ranged from 0.44 to 0.57. Exact and within-category agreement ranged from 41.4%-50.8% and 72.0%-79.0%. Baseline-referenced LWK ranged from 0.22 to 0.47, with Spearman correlations of 0.29 to 0.56. RPAF correlated with total IPAQ leisure score (rs = 0.60), IPAQ walking score (rs = 0.58), SF-36-PF (rs = 0.33) and TUG score (rs = -0.25). No significant correlation was identified between RPAF, HGS or WhTR. RPAF discriminated known groups for WHO guideline-sufficient activity, SF-36-PF, and TUG fall risk. The RPAF item demonstrated fair-to-moderate agreement in consecutive waves, with weaker baseline-referenced tracking. Cross-sectional validity was highest with total IPAQ leisure activity. The item may provide a pragmatic measure for RPAF in womens cohort studies.

Background: Aerobic exercise is recommended during stroke rehabilitation to improve cardiorespiratory fitness and support recovery; however, participation rates remain low. While institutional and system-level barriers have been widely examined, less is known about how individual patient factors influence engagement in aerobic exercise during rehabilitation. Objectives: We aimed to determine whether depressive symptoms, apathy, self-efficacy and outcome expectations for exercise, perceived barriers, or past exercise history were associated with aerobic exercise participation in stroke rehabilitation. Methods: In this prospective cohort sub-study, adults admitted to in- or out-patient stroke rehabilitation at three urban hospitals completed validated questionnaires assessing depressive symptoms, apathy, exercise self-efficacy, outcome expectations for exercise, perceived barriers to being active, and premorbid exercise history. Participants were separated into two groups for analysis: those who completed aerobic exercise during rehabilitation and those who did not. Equivalence testing and between-group comparisons were performed. Results: Sixty-two participants were enrolled; 16 participated in aerobic exercise and 46 did not. Groups were not equivalent on any individual-level factors. Compared to non-participants, those who performed aerobic exercise had significantly higher depressive symptom scores (p=0.0025) and lower self-efficacy for exercise (p=0.0087). Non-participants demonstrated significantly higher apathy (p=0.0007). No significant differences were found for outcome expectations, perceived barriers, or exercise history. Conclusion: Depressive symptoms and lower self-efficacy did not impede aerobic exercise participation during rehabilitation. Increased apathy, however, was associated with non-participation. Findings highlight the need for individually tailored aerobic exercise prescriptions that consider motivational and affective factors to optimize engagement during stroke rehabilitation.

Abstract Preventable neonatal deaths remain a major public health problem in Zimbabwe, where near-universal antenatal and facility-delivery coverage coexist with a rising neonatal mortality rate. This study evaluated whether institutionalising three core "vital signs" of the community health system (a trained village health worker (VHW) workforce, functional community governance structures, and modified women's and men's participatory learning and action groups) reduces preventable neonatal deaths in Mashonaland West Province. An embedded QUAN (qual) mixed-methods design was used, with a two-arm, parallel-group cluster-randomised controlled trial as the dominant strand. Fifty-two ward-level clusters were randomised 1:1 to the institutionalised community health system package or to standard Ministry of Health and Child Care community services, and 984 pregnant women were enrolled between 1 September 2020 and 31 October 2021, with each mother-infant pair followed to 28 days after delivery, yielding 973 mother-infant pairs for intention-to-treat analysis. The primary outcome was neonatal death within 28 days of life, expressed per 1,000 live births. The primary analysis used a three-level mixed-effects log-binomial regression model with cluster and community-health-worker random intercepts, adjusted for pre-specified covariates. Supervised machine-learning classifiers with leave-one-cluster-out cross-validation, Cox proportional-hazards regression, and multilevel logistic models were fitted as supplementary analyses. An embedded longitudinal process evaluation used key informant interviews and focus group discussions, which were analysed thematically and integrated with the quantitative findings. The neonatal mortality rate was 44.8 per 1,000 live births in the intervention arm versus 110.1 per 1,000 in the control arm. The adjusted risk ratio for neonatal death was 0.43 (95% CI 0.26-0.70; p < 0.001), a 57% relative reduction, with a number needed to treat of 16 mother-infant pairs (95% CI 11-29). Low birthweight (<2,500 g), birth interval under two years, and low community women's literacy were the strongest risk factors, while trained VHWs, functional community governance, early antenatal care, and sustained participatory group attendance were independently protective. The women's and men's groups were protective in a dose-dependent manner, becoming significant at four or more cycles (about 14 meetings) (adjusted odds ratio 0.71; 95% CI 0.60-0.85; p = 0.001). A random forest classifier discriminated against neonatal deaths with a cross-validated area under the curve of 0.842 and a sensitivity of 0.912. Qualitative findings converged with the trial results, identifying male engagement, earlier care-seeking, danger-sign literacy, social-network activation, and community death audits as the behavioural and structural mechanisms of change. Institutionalising the community health system package (trained VHWs, functional governance, early antenatal engagement, and sustained participatory groups) was associated with a substantial reduction in preventable neonatal deaths. The findings suggest that in high-coverage, high-mortality settings, the binding constraint is structural rather than clinical, and that scaling functional community governance and workforce infrastructure in the most disadvantaged communities may accelerate progress toward neonatal survival targets. The principal limitations are a one-year follow-up period, the rarity of neonatal death, and concurrent national programming that only partially reached the control clusters. Trial registration: Pan African Clinical Trials Registry, PACTR202607591142118 (https://pactr.samrc.ac.za/TrialDisplay.aspx?TrialID=PACTR202607591142118); registered retrospectively on 7 July 2026.

Biological aging occurs heterogeneously across individuals and organs. However, current measures of biological age incompletely capture organ-specific differences in health and disease risk. Because chest CT visualizes multiple thoracic organs, it offers an opportunity to quantify structural aging across organ systems. Here, we developed MOSAIC-Age, a framework characterizing eight organ-specific aging clocks on chest CT. The clocks were developed and validated using 9,971 CT scans from CT-RATE and MIDRC, and subsequently locked and applied to two independent prospective cohorts with 35,293 participants from the National Lung Screening Trial and Genetic Epidemiology of COPD study. CT-derived biological age gaps (BAGs) were examined in relation to lifestyle and socioeconomic factors, prevalent comorbidities, incident chronic diseases, and all-cause and cause-specific mortality. Higher BAGs, indicating organs that appeared older on CT than expected for their chronological age, were broadly associated with adverse health characteristics, chronic disease burden, and increased mortality risk. Multiple disease outcomes were associated with aging across several organs, whereas in multivariable analyses including all eight organ-specific BAGs, the remaining associations were more organ specific. A greater number of markedly older-appearing organs and a faster pace of aging were each associated with higher mortality. Together, these findings demonstrate that routine chest CT captures both shared and organ-specific patterns of biological aging and establish CT-derived organ aging as a quantitative imaging biomarker for assessing multi-organ health and long-term disease risk.

Background Rural-urban disparities in SARS-CoV-2 infection epidemiology remain poorly quantified and understood in Zambia despite differences in healthcare access, services and preventive interventions. This study examined the geographical distribution and associated factors of SARS-CoV-2 cases across selected rural and urban districts of Zambia. Methods A convergent mixed-methods study comprised of quantitative survey and qualitative interviews was conducted in; Ndola (Urban), Kafue (Peri-urban) and Lufwanyama (Rural). The proximate determinant framework guided variable selection and interpretation. Quantitative combined (Hospital-surveillance data with community survey), while qualitative included In-depth interviews. Participants were sampled using multistage sampling technique. Quantitative data were analysed using STATA version 17, while qualitative data were analysed thematically. Findings were integrated through triangulation. Results A total of 528 participants were included, with a median age 31 years (15-71). Overall SARS-CoV-2 positivity was 12.6%, varying across rural (16.5%), peri-urban (14.9%), and urban (9.9%) settings, though residence was not associated with infection (P<0.132). Participants aged [&ge;]49 years had significantly higher odds of infection (aOR=8.78; 95% CI:1.15-66.99), whereas secondary education (aOR=0.37; 95% CI:0.16-0.86) and hospital-based testing (aOR=0.37; 95% CI:0.15-0.92) were associated with lower odds of infection. Vaccine uptake was highest in urban areas but was not independently associated with infection. Qualitative findings revealed marked rural-urban differences in perceived susceptibility, testing access, vaccine decision-making, and adherence to preventive measures, explaining several quantitative observations. Conclusion SARS-CoV-2 infection across rural and urban settings in Zambia was influenced by demographic, behavioral, and health-system factors rather than geographic residence alone. These findings highlight the need for context-specific prevention strategies, equitable access to testing, strengthened community surveillance, and targeted risk communication to improve preparedness and response for future respiratory disease outbreaks.